CFTP/12-004 
PCCF RI 12-03 



Revisiting the T (K — > eu) /Y (K — > fiu) ratio in 
supersymmetric unified models 

Renato M. Fonseca", J. C. Romao" and A. M. Teixeira 6 



a Centro de Ffsica Teorica de Parti'culas, CFTP, Instituto Superior Tecnico, 
Universidade Tecnica de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 

b Laboratoire de Physique Corpusculaire, CNRS/IN2P3 - UMR 6533, 
Campus des Cezeaux, 24 Av. des Landais, F-63171 Aubiere Cedex, France 



Abstract 

It has been pointed out that supersymmetric extensions of the Standard Model can induce 
significant changes to the theoretical prediction of the ratio T (K — > ev) /T (K — > jiv) = Rk, 
through lepton flavour violating couplings. In this work we carry out a full computation of all 
one-loop corrections to the relevant u£H + vertex, and discuss the new contributions to Rk 
arising in the context of different constrained (minimal supergravity inspired) models which 
succeed in accounting for neutrino data, further considering the possibility of accommodating 
a near future observation of a \i — > ey transition. We also re-evaluate the prospects for 
Rk in the framework of unconstrained supersymmetric models. In all cases, we address the 
question of whether it is possible to saturate the current experimental sensitivity on Rk while in 
agreement with the recent limits on B-meson decay observables (in particular BK(B S — > 
and BR(£>„ — > tv)), as well as BR(r — > ey) and available collider constraints. Our findings 
reveal that in view of the recent bounds, and even when enhanced by effective sources of 
flavour violation in the right-handed e — f sector, constrained supersymmetric (seesaw) models 
typically provide excessively small contributions to Rk- Larger contributions can be found 
in more general settings, where the charged Higgs mass can be effectively lowered, and even 
further enhanced in the unconstrained MSSM. However, our analysis clearly shows that even in 
this last case SUSY contributions to Rk are still unable to saturate the current experimental 
bounds on this observable, especially due to a strong tension with the B u — > tv bound. 
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1 Introduction 



Neutrino oscillations have provided the first experimental manifestation of flavour violation in the 
lepton sector, fuelling the need to consider extensions of the Standard Model (SM) that succeed 
in explaining the smallness of neutrino masses and the observed pattern of their mixings [l|- 
[3]. In addition to the many facilities dedicated to study neutral leptons, there is currently a 
great experimental effort to search for signals of flavour violation in the charged lepton sector 
(cLFV), since such an observation would provide clear evidence for the existence of new physics 
beyond the SM (trivially extended to accommodate massive neutrinos). The quest for the origin 
of the underlying mechanism of flavour violation in the lepton sector has been actively pursued in 
recent years, becoming even more challenging as the MEG experiment is continually improving the 
sensitivity to fi — > decays |4|, thus opening the door for a possible measurement (observation) in 
the very near future. The current bounds on other radiative decays (i.e. l{ — > ^7), or three-body 
decays [li — > Mj) are already impressive |5], and are expected to be further improved in the future. 

Supersymmetric (SUSY) extensions of the SM offer new sources of CP and flavour violation, in 
both quark and lepton sectors. Given the strong experimental constraints, especially on CP and 
flavour violating observables involving the strongly interacting sector, phenomenological analyses 
in general favour the so-called "flavour-blind" mechanisms of SUSY breaking, where universality 
of the soft breaking terms is assumed at some high energy scale: in these constrained scenarios, 
the only sources of flavour violation (FV) are the quark and charged lepton Yukawa couplings. 
In order to accommodate current neutrino data, mechanisms of neutrino mass generation, such 
as the seesaw (in its different realisations - for a review of the latter, see for instance (6][tJ ) , are 
often implemented in the framework of (constrained) SUSY models: in the case of the so-called 
"SUSY-seesaw", radiatively induced flavour violation in the slepton sector (8) can provide sizable 
contributions to cLFV observables. The latter have been extensively studied, both at high- and 
low-energies, over the past years (see e.g. [9]). Flavour violation can be also incorporated in a 
more phenomenological approach, where at low-energies new sources of FV are present in the soft 
SUSY breaking terms. However, these are severely constrained by a large number of observables 



(see, e.g. 10 and references therein). 

In addition to the above mentioned rare lepton decays, leptonic and semi-leptonic meson decays 
also offer a rich testing ground for cLFV. Here we will be particularly interested in leptonic K 
decays, which (as is also the case of leptonic 7r decays) constitute very good probes of violation of 
lepton universality. The potential of these observables, especially regarding SUSY extensions of 



the SM, was firstly noticed in and later investigated in greater detail in [l2 - 14 



By themselves, these decays are heavily hampered by hadronic uncertainties and, in order to 
reduce the latter (and render these decays an efficient probe of new physics), one usually considers 
the ratio 

R * = r(K+^M7])' ( 5 

since in this case the hadronic uncertainties cancel to a very good approximation. As a consequence, 



the SM prediction can be computed with high precision 15-17 . The most recent analysis has 
provided the following value (17): 



Rj^ 1 = (2.477 ± 0.001) x 10" 5 . (1.2) 
On the experimental side, the NA62 collaboration has recently obtained very stringent bounds fl8] : 

R e ^ = (2.488 ±0.010) x 10~ 5 , (1.3) 
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which should be compared with the SM prediction (Eq. (1.2)). In order to do so, it is often useful 
to introduce the following parametrisation, 

R e *P = Rf 1 (1 + Ar) , Ar = Rk/bP* - 1 , (1.4) 

where Ar is a quantity denoting potential contributions arising from scenarios of new physics (NP). 



Comparing the theoretical SM prediction to the current bounds (i.e., Eqs. (1.2 1.3)), one verifies 
that observation is compatible with the SM (at la) for 

Ar = (4 ±4) x 1(T 3 . (1.5) 

Previous analyses have investigated supersymmetric contributions to Rk in different frame- 
works, as for instance low-energy SUSY extensions of the SM (i.e. the unconstrained Minimal Su- 
persymmetric Standard Model (MSSM)) [IT 12 14j, or non-minimal grand unified models (where 



higher dimensional terms contribute to fermion masses) |13| . These studies have also considered 
the interplay of Rk with other important low-energy flavour observables, magnetic and electric 
lepton moments and potential implications for leptonic CP violation. Distinct computations, based 
on an approximate parametrisation of flavour violating effects - the Mass Insertion Approximation 
(MIA) |19| - allowed to establish that SUSY LFV contributions can induce large contributions to 
the breaking of lepton universality, as parametrised by Ar. The dominant FV contributions are in 
general associated to charged-Higgs mediated processes, being enhanced due to non-holomorphic 



effects - the so-called "HRS" mechanism 20 -, and require flavour violation in the RR block of 
the charged slepton mass matrix. It is important to notice that these Higgs contributions have 
been known to have an impact on numerous observables, and can become especially relevant for 



the large tan /3 regime 20 - 31 



In the present work, we re-evaluate the potential of a broad class of supersymmetric extensions 
of the SM to saturate the current measurement of Rk- Contrary to previous studies, we conduct 
a full computation of the one-loop corrections to the vlH + vertex, taking into account the im- 
portant contributions from non-holormophic effective Higgs-mediated interactions. When possible 
we establish a bridge between our results and approximate analytical expressions in the litera- 
ture, and we stress the potential enhancements to the total SUSY contributions. In our numerical 
analysis we re-investigate the prospects regarding Rk of a constrained MSSM onto which several 
seesaw realisations are embedded (type I |32) and II |33] , as well as the inverse seesaw |34] ) , also 



briefly addressing L-R symmetric models 35 36 . We then consider more relaxed scenarios, such 



as non-universal Higgs mass (NUHM) models at high-scale (which are known to enhance this class 
of observables |13| due to potentially lighter charged Higgs boson masses) , and discuss the general 
prospects of unconstrained low-energy SUSY models. In all cases, we revisit the Rk observable 
in the light of new experimental data: in addition to LHC bounds^] on the sparticle spectrum 



and a number of low-energy flavour-related bounds |4l 5l, we implement the very recent LHCb 



38 



results concerning the BR(S S — >• /j + fi~) 39 . As we discuss here, the increasing tension with 
low-energy observables, in particular with B u —> tv, precludes sizable SUSY contributions to Rk 
even in the context of otherwise favoured candidate models as is the case of semi-constrained and 
unconstrained SUSY models. 

This document is organised as follows. Section [2] is devoted to the computation of the 1-loop 
MSSM prediction for Rk- We compare our (full) result to the approximations in the literature 



In our numerical analysis we do not require the lightest Higgs to be in strict agreement with recent LHC search 
results |37) : while in the general the case (especially for constrained (seesaw) models), we only favour regimes where 
its mass is larger than 118 GeV, when considering semi-constrained and unconstrained models, a significant part of 
the studied region does indeed comply with rrih ~ 125 GeV. 
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by means of the mass insertion approximation (among other simplifications), and discuss the 
dominant sources of flavour violation, and the implications to other observables. Our results for a 
number of models are collected in Section [3] Further discussion and concluding remarks are given 
in Section [4] In the Appendices, we detail the computation of the renormalised charged lepton - 
neutrino - charged Higgs vertex, and summarise the key features of two supersymmetric seesaw 
realisations (types I and II) used in the numerical analysis. 



2 Supersymmetric contributions to Rk 

In the SM, the decay widths of pseudoscalar mesons into light leptons are given by 



G 2 F m P mj 
8vr 



1 



m, 



fpWr. 



gq 



(2.1) 



where P denotes n,K,D or B mesons, with mass mp and decay constant fp, and where Gf is 
the Fermi constant, me the lepton mass and V qq i the corresponding Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element. These decays are helicity suppressed (as can be seen from the factor mf in 



Eq. (2.1)), and the prediction for their amplitude is thus hampered by the hadronic uncertainties 
in the meson decay constants. As mentioned in the Introduction, ratios of these amplitudes are 
independent of fp to a very good approximation, and the SM prediction can then be computed 
very precisely. Concerning the kaon decay ratio Rk, the SM prediction (inclusive of internal 
bremsstrahlung radiation) is |17] 



R K 




mz 



mf 



(1 + 5R QEB ) 



(2.2) 



where 5Rqed is a small electromagnetic correction accounting for internal bremsstrahlung and 
structure-dependent effects (5R QE d = (-3.60 ± 0.04)% 





H+ 



K 




Figure 1: Tree level contributions to Rk - SM and charged Higgs. 
In supersymmetric models, the extended Higgs sector can play an important role in lepton 



flavour violating transitions and decays (see 20 -31]). The effects of the additional Higgs are also 
sizable in meson decays through a charged Higgs boson, as schematically depicted in Fig. [T] In 
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particular, for kaons, one finds |21] 



l- t an 2 /3^ ms ) : (2.:!) 
m s + m u J 

however, despite this new tree-level contribution, Rk is unaffected, as the extra factor does not 
depend on the (flavoured) leptonic part of the process. 

New contributions to Rk only emerge at higher order: at one-loop level, there are box and 
vertex contributions, wave function renormalisation, which can be both lepton flavour conserv- 
ing (LFC) and lepton flavour violating. Flavour conserving contributions arise from loop cor- 
rections to the propagator, through heavy Higgs exchange (neutral or charged) as well as 
from chargino/neutralino-sleptons (in the latter case stemming from non- universal slepton masses, 



i.e., a selectron-smuon mass splitting). As concluded in 11 , in the framework of SUSY models 
where lepton flavour is conserved, the new contributions to Ar SUSY are too small to be within 
experimental reach. 

On the other hand, Higgs mediated LFV processes are capable of providing an important 
contribution when the kaon decays into a electron plus a tau-neutrino. For such LFV Higgs 
couplings to arise, the leptonic doublet (L) must couple to more than one Higgs doublet. However, 
at tree level in the MSSM, L can only couple to Hi, and therefore such LFV Higgs couplings arise 
only at loop level, due to the generation of an effective non-holomorphic coupling between L and 
H2 - the HRS mechanism |20| - which is a crucial ingredient in enhancing the Higgs contributions 
to LFV observables. In what follows, we address the impact of these non-holomorphic terms for 
Rk- 

2.1 LFV Higgs mediated contributions to Rk 

We consider as starting point the MSSM, defined by its superpotential and soft-SUSY breaking 
Lagrangian. We detail below the relevant terms for our discussion: 

W = U C Y U QH 2 - D c Y d QHi - E c Y l LH\ - ^HiH 2 , (2.4) 

Vsoft = - Aloft = ( M a 4>a + h..C.) + TTlg. H* Hi 

+ (B Hi H 2 + h.c.) + I* L m\ i L + t R m\l R 

+ {H x t R A l i L +h.c.) + ..., (2.5) 

where M a denotes the soft-gaugino mass terms, stand for the squark terms, and we have 

omitted flavour indices. For the SU(2) superfield products, we adopt the convention Hi H 2 = 
H\ H\ — Hf H\ (and likewise for similar cases). 

From an effective theory approach, the HRS mechanism can be accounted for by additional 
terms, corresponding to the higher-order corrections to the Higgs-neutrino-charged lepton interac- 
tion (schematically depicted in Fig. [2| . 

At tree-level, the Lagrangian describing the vtH^ 1 interaction is given by 

Cfi* =V L Y l U R Hi* + h.c. 

= (2 3 / 4 G^ /2 ) tan V L M l t R H + + h.c. , (2.6) 
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Figure 2: Corrections to the v£H + vertex, as discussed in the text. 

with M l = diag (m e , m„, m T ). At loop level, two new terms are generated: VlA^ IrH^ — 
JlA°£jiH2 + h.c. The second one, with A , forces a redefinition of the charged lepton Yukawa 

couplings, y't = if > Y 1 ^ rs ^ A°tan/3, which in turn implies a redefinition of the charged 

lepton propagator; the term with A + corrects the Higgs-neutrino-charged lepton vertex^] Once 



these terms are taken into account, the interaction Lagrangian, Eq. (2.6), becomes 

C H± = (2 3 / 4 G 1 ' 2 ) t a n(3v L M l e R H + 



+ cos p V L ( A+ - A tan 2 0) £ R H + + h.c. . (2.7) 

Since in the SU(2) L -preserving limit A+ = A , it is reasonable to assume that, after electroweak 
(EW) symmetry breaking, both terms remain approximately of the same order of magnitude. 
Hence, it is clear that the contribution associated with A (the loop contribution to the charged 
lepton mass term) will be enhanced by a a factor of tan 2 j3 when compared to the one associ- 
ated with A + . This simple discussion allows to understand the origin of the dominant SUSY 
contribution^] to Rr. 

As we proceed to discuss, a quantitative assessment of the corrections to A + and A requires 
considering the higher-order effects on the vertex Vl Z h £rH + (see also |4l|). The Z H matrix 
depends on the following (loop-induced) quantities: 

• rf L and r\ u L (corrections to the kinetic terms of £l and vl)\ 

• rf m (correction to the charged lepton mass term); 

• r\ H (correction to the u£H vertex). 

The expressions for the distinct ^-parameters can be found in Appendix [a| Instead of Z H , which 
includes both tree and loop level effects, it proves to be more convenient to use the following 
combination, 

tan^ /m, V ^n^_ ^ / a-i ^ £l + 
L l 2 \m H + J m s + m u V / 



23/4G F 



where 



-tan 2 /3 ^L) ms , (2.9) 



A 




m H + J m s +m u 



2 3 / 4 G X l 2 tan/3 




(2.10) 
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An extensive discussion on the radiatively induced couplings which are at the origin of the HRS effect can be 



found in 



40 



3 There are additional corrections to the qq'H ± vertex, which are mainly due to a similar modification of the the 
quark Yukawa couplings - especially that of the strange quarks. This amounts to a small multiplicative effect on 
Ar which we will not discuss here (see [14[ for details). 
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In the above, e encodes the tree level Higgs mediated amplitude (which does not change the SM 
prediction for Rk), while A, a matrix in lepton flavour space, encodes the 1-loop effects. The main 
contribution is expected to arise from ff m . 

The Ar observable is then related to e and A as follows: 



Ar 



Rk 
Rf 1 



1 + 



At 
1+e 



1 + 



1+e 



1 + T 



At 



A 

1+e 



1. 



If the slepton mixing is sufficiently large, this expression can be approximated as 

Ar « 2 Re(A ee ) + (A t A) ee . 



(2.11) 



(2.12) 



In the above, the first (linear) term on the right hand-side is due to an interference with the SM 
process, and is thus lepton flavour conserving. As shown in pT| , this contribution can be enhanced 
through both large RR and LL slepton mixing. On the other hand, the quadratic term (A^A) ee 
can be augmented mainly through a large LFV contribution from A re , which can only be obtained 
in the presence of significant RR slepton mixing. 



2.2 Generating Ar: sources of flavour violation and experimental constraints 

In order to understand the dependence of Ar on the SUSY parameters, and the origin of the 
dominant contributions to this observable, an approximate expression for A is required. Firstly, we 



notice that the previous discussion, leading to Eq. (2.7), suggests that the rf m term is responsible for 



the dominant contributions to Ar. Thus, in what follows, and for the purpose of obtaining simple 
analytical expressions, we shall neglect the contributions of the other terms (although these are 
included in the numerical analysis of Section [3J . A fairly simple analytical insight can be obtained 
when working in the limit in which the virtual particles in the loops (sleptons and gauginos) are 
assumed to have similar masses, so that their relative mass splittings are indeed small. In this 
limit, one can Taylor-expand the loop functions entering rf m (see Appendix Ar working to third 



order in this expansion, and keeping only the terms enhanced by a factor of m T tan (3 (where 
^SUSY) ttiew denote the SUSY breaking scale and EW scale, respectively), we obtain 



Ar 



l + Xll 



10 ml „ 




-fj, 2 + S 3 



3 fi 2 + 2Ml 
10 



— 2 



(2.13) 



where /i, Mi and (m~) eT denote the low-energy values of the Higgs bilinear term, bino soft-breaking 
mass, and off-diagonal entry of the soft-breaking left-handed slepton mass matrix, respectively. We 
have also introduced m~ = \ ((m|) + (m 2 )), the average mass squared of sleptons and neutralinos 



m SUSY)' and 5 



X 



(m 2 )), the corresponding splitting. The quantity X is given by 

X 



1927T 



2/2 I ^ tan 3 f3 m 
2 m K g uMi 



r H). 



in 



H+ 



(m L° )3 



(2.14) 



and it illustrates in a transparent (albeit approximate) way the origin of the terms contributing 
to the enhancement of Rk ■ in addition to the factor tan 3 /5 / m? H+ , usually associated with Higgs 
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exchanges, the crucial flavour violating source emerges from the off-diagonal (re) entry of the 
right-handed slepton soft-breaking mass matrix. 

Using the above analytical approximation, one easily recovers the results in the literature, 
usually obtained using the MIA. For instance, Eq. (11) of Ref. Ill] amounts to 



Ar ~ 2X (ml) + X 2 (mfj 2 + X 2 5 2 , (2.15) 



which stems from having kept the dominant (crucial) second and third order contributions in the 
expansion: X 2 5 2 and 2X I mf ) + X 2 I m~ ) , respectively. 



Regardless of the approximation considered, it is thus clear that the LFV effects on kaon decays 
into a ev or [iv pair can be enhanced in the large tan j3 regime (especially in the presence of low 

values of m H +), and via a large RR slepton mixing ( m ^J ■ Although the latter is indeed the 



privileged source, notice that, as can be seen from Eq. (2.15 ), a strong enhancement can be obtained 



from sizable flavour violating entries of the left-handed slepton soft-breaking mass, y m L j ■ This 
is in fact a globally flavour conserving effect (which can also account for negative contributions to 
Rk)- Previous experimental measurements of Rk appeared to favour values smaller than the SM 



theoretical estimation, thus motivating the study of regimes leading to negative values of Ar 11 



but these regimes have now become disfavoured in view of the present bounds, Eq. (1.5). 

Clearly, these Higgs mediated exchanges, as well as the FV terms at the origin of the strong 
enhancement to Rk, will have an impact on a number of other low-energy observables, as can be 



easily inferred from the structure of Eqs. ( 2,13{[2~T5 ). This has been extensively addressed in the 



literature |ll]|l4| , and here we will only briefly discuss the most relevant observables: electroweak 
precision data on the anomalous electric and magnetic moments of the electron, as well as the 
naturalness of the electron mass, directly constrain the rf m corrections (and rf L , rj H ); low-energy 
cLFV observables, such as r — > and r — > 31 decays are also extremely sensitive probes of Higgs 
mediated exchanges, and in the case of r — e transitions, depend on the same flavour violating 
entries. It has been suggested that positive and negative values of Ar can be of the order of 



1%, still in agreement with data on the electron's electric dipole moment and on r — > I'y 11-13 
Finally, other meson decays, such as B — > li (and B — > iv), exhibit a similar dependence on tan/3, 
tan"/3/m^+ (42] (n ranging from 2 to 6, depending on the other SUSY parameters), and may 
also lead to indirect bounds on Ar. In particular, the strict bounds on BR(S U — > tu) [5] and 



the very recent limits on BR(i? s — > fi + fi ) 39 might severely constrain the allowed regions in 



SUSY parameter space for large tan /3. Although we will come to this issue in greater detail when 
discussing the numerical results, it is clear that the similar nature of the K + — > Iv and B u — > tv 



processes (easily inferred from a generalization of Eq. (2.3), see e.g. |2H|43| ) will lead to a tension 
when light charged Higgs masses are considered to saturate the bounds on Rk- 

Supersymmetric models of neutrino mass generation (such as the SUSY seesaw) naturally 
induce sizable cLFV contributions, via radiatively generated off-diagonal terms in the LL (and to 
a lesser extent LR) slepton soft-breaking mass matrices p]. In addition to explaining neutrino 
masses and mixings, such models can also easily account for values of BR(/x — > e^f), within the 
reach of the MEG experiment. In view of the recent confirmation of a large value for the Chooz 
angle (#13 ~ 8.8°) |3| and on the impact it might have on (m|) er , in the numerical analysis of the 



following section we will also consider different realisations of the SUSY seesaw (type I [32] , II 33 
and inverse |34| ) , embedded in the framework of constrained SUSY models. We will also revisit 
semi-constrained scenarios allowing for light values of mjj+, re-evaluating the predictions for Rk 
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under a full, one loop-computation, and in view of recent experimental data. Finally, we confront 
these (semi-)constrained scenarios with general, low-energy realisations, of the MSSM. 



3 Prospects for Rx' unified vs unconstrained SUSY models 

In this section we evaluate the SUSY contributions to Rk, with the results obtained via the full 
expressions for Ar, as described in Section [2] These were implemented into the SPheno public 
code |44) , which was accordingly modified to allow the different studies. It is important to stress 
that although some approximations have still been done (as previously discussed) , the results based 
on the present computation strongly improve upon those so far reported in the literature (mostly 
obtained using the MIA). Although the different contributions cannot be easily disentangled due 
to having carried a full computation, our results automatically include all one-loop lepton flavour 
violating and lepton flavour conserving contributions (in association with charged Higgs mediation, 
see footnote [3]). As mentioned before, we evaluate Rk in the framework of constrained, semi- 
constrained (NUHM) and unconstrained SUSY models. Concerning the first two, we assume 
some flavour blind mechanism of SUSY breaking (for instance minimal supergravity (mSUGRA) 
inspired), so that the soft breaking parameters obey universality conditions at some high-energy 
scale, which we choose to be the gauge coupling unification scale Mqut ~ 10 16 GeV, 



{m D f.. = (m D ) 2 .. = (m L )l = (m A )* 



: m 



o v ij i 



A 1 ) =A (Y%. (3.1) 

/ ij 

In the above, tuq and Aq are the universal scalar soft-breaking mass and trilinear couplings of 
the cMSSM, and i,j denote lepton flavour indices (z,j = 1,2,3). In the latter case, the gaugino 
masses are also assumed to be universal, their common value being denoted by Mi/ 2 . We will 
also consider the supersymmetrisation of several mechanisms for neutrino mass generation. More 
specifically, we have considered the type I and type II SUSY seesaw (as detailed in Appendix [B|) . 
We briefly comment on the inverse SUSY seesaw, and discuss a L — R model. 



The strict universality boundary conditions of Eqs. (3.1) will be relaxed for the Higgs sector 



when we address NUHM scenarios, so that in the latter case we will have 

m 2 Hl ± m\ 2 j= ml . (3.2) 

All the above universality hypothesis will be further relaxed when, for completeness, and to allow 
a final comparison with previous analyses, we address the low-energy unconstrained MSSM. 



In our numerical analysis, we took into account LHC bounds on the SUSY spectrum |38j, as 
well as the constraints from low-energy flavour dedicated experiments 
In particular, concerning lepton flavour violation, we have considered 



, and neutrino data [lj|2j. 
151: 



BR(r -> ery) < 3.3 x 1(T 8 (90%C.L.) , (3.3) 

BR(r -> 3 e) < 2.7 x 1(T 8 (90%C.L.) , (3.4) 

BR(/j -»■ &y) < 2.4 x 1(T 12 (90%C.L.) , (3.5) 

BR(B U -> tu) > 9.7 x 1(T 5 (2 a). (3.6) 
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m (GeV) 


Mi/a (GeV) 


tan j3 


A) (GeV) 


sign(/x) 


10.3.1 


300 


450 


10 





1 


P20 


330 


500 


20 


-500 


1 


P30 


330 


500 


30 


-500 


1 


40.1.1 


330 


500 


40 


-500 


1 


40.3.1 


1000 


350 


40 


-500 


1 



Table 1: cMSSM (benchmark) points used in the numerical analysis. 
Also relevant are the recent LHCb bounds 1 39 1 



BR{B S ->• fi + fi-) < 4.5 x 10~ 9 (95%GL.), (3.7) 
BR{B -»• n+fjT) < 1.03 x 10~ 9 (95%C.L.) . (3.8) 

When addressing models for neutrino mass generation, we take the following (best-fit) values 
for the neutrino mixing angles [2] (where #13 is already in good agreement with the most recent 
results from |3]), 

sin 2 12 = 0.312±°;°& sin 2 23 = 0.52±g;gf, 

S in 2 13 ~O.O13i°;°oL ( 3 - 9 ) 
Am 2 2 = (7.59i°; 2 °) x 10~ 5 eV 2 , (3.10) 

Am 2 3 = (2.50±g:5§) x 10_3 ey2 - ( 3 - n ) 

Regarding the leptonic mixing matrix (£/mns) we adopt the standard parametrisation. In the 
present analysis, all CP violating phases are set to zercQ 

3.1 mSUGRA inspired scenarios: cMSSM and the SUSY seesaw 

We begin by re-evaluating, through a full computation of the one-loop corrections, the maximal 
amount of supersymmetric contributions to Rk in constrained SUSY scenarios. For a first evalua- 
tion of Rk, we consider different cMSSM (mSUGRA-like) points, defined in Table [l] Among them 
are several cMSSM benchmark points from [45 j, representative of low and large tan/3 regimes, 
as well as some variations. Notice that, as mentioned before, these choices are compatible with 
having a Higgs boson mass above 118 GeV but will be excluded once we require rrih to lie close to 
125 GeV as suggested by LHC results [37]. 



As could be expected from Eqs. (2.13 2.15), in a strict cMSSM scenario (in agreement with 
the experimental bounds above referred to) the SUSY contributions to Rk are extremely small; 
motivated by the need to accommodate neutrino data, and at the same time accounting for values 
of BR(/j —> cy) within MEG reach, we implement type I and type II seesaws in mSUGRA- 



inspired models (see Appendices B.l and B.2 ). Regarding the heavy-scale mediators, we considered 
degenerate right-handed neutrinos, as well as degenerate scalar triplets. We set the seesaw scale 
aiming at maximising the (low-energy) non-diagonal entries of the soft-breaking slepton mass 



matrices, while still in agreement with the current low-energy bounds (see Eqs. (3.3 3.8)). In 



particular, we have tried to maximise the LL contributions to Ar, i.e., (m~) eT , and to obtain 



4 We will assume that we are in a strictly CP conserving framework, where all terms are taken to be real. This 
implies that there will be no contributions to observables such as electric dipole moments, or CP asymmetries. 
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BR(/i ->■ cy) within MEG reach (i.e. 10~ 13 < BR(/j -> e7)< 2.4 x 10~ 12 ). However, and due 
to the fact that both seesaw realisations fail to account for radiatively induced LFV in the right- 
handed slepton sector, one finds values |Ar| < 2 x 10 -8 . It is worth emphasising that if one further 
requires nih to lie close to 125 GeV (as suggested by recent findings |37|), then one is led to regions 
in mSUGRA parameter space where, due to the much heavier sparticle masses and typically lower 
values of tan/3, the SUSY contributions to Rx become even further suppressed. 

Thus, and even under a full computation of the corrections to the vlH + vertex, we nevertheless 



confirm that, as firstly put forward in the analyses of 11 12 strictly constrained SUSY and SUSY 



seesaw models indeed fail to account for values of Rk close to the present limits. 

Clearly, new sources of flavour violation, associated to the right-handed sector are required: in 

what follows, we maintain universality of soft-breaking terms allowing, at the grand unified (GUT) 

scale, for a single r — e flavour violating entry in m\. This approach is somewhat closer to the 
i — 1 1 — i 

lines of 111 - 14 , although in our computation we will still conduct a full evaluation of the distinct 
contributions to Ar, and we consider otherwise universal soft-breaking terms. Without invoking a 
specific framework/scenario of SUSY breaking that would account for such a pattern, we thus set 

6«« = ^^0. (3.12) 

As discussed above, low-energy constraints on LFV observables (especially r — > cy), severely 
constrain this entry. 

In Fig. |3j we present our results for Ar scanning the mo — Afi/2 plane for a regime of large 
tan/3. We have set 8§f = 0.1, tan/3 = 40, and taken A = -500 GeV. The surveys displayed 
in the panels correspond to having embedded a type I (left) or type II (right) seesaw onto this 
near-mSUGRA framework. 

As can be readily seen from Fig. |3j once the constraints from low-energy observables have been 
applied, in the type I SUSY seesaw, the maximum values for Ar are 0(1O -7 ), associated to the 
region with a lighter SUSY spectra (which is in turn disfavoured by a "heavy" light Higgs). Even for 
the comparatively small non- universality, S^f* = 0.1, a considerable region of the parameter space 
is excluded due to excessive contributions to BR(B U — > tv) and BR(r — > ej), thus precluding the 
possibility of large values of Ar. In a regime of large tan/3, the contributions to BR(5 S — > fj> + /j>~) 
are also sizable, and the recent LHCb results seem to exclude the regions of the parameter space 
where one could still have Ar ~ O(10" 6 '" 7 ). The excessive SUSY contributions to BR(i? s — > 
can be somewhat reduced by adjusting Aq (in Fig. [3] we fixed Aq = —500 GeV) and the 
values of Ar can be slightly augmented by increasing S^; in the latter case, the r — > bound 
proves to be the most constraining, and values of Ar larger than 0(1O~ 6 ' -7 ) cannot be obtained 
in these constrained SUSY seesaw models. 

The situation is somewhat different for the type II case: firstly notice that a sizable region in 
the mo — M1/2 plane is associated to negative contributions to Rk, which are currently disfavoured. 
In the remaining (allowed) parameter space, the values of Ar are slightly smaller than for the type 
I case: this is a consequence of a non trivial interplay between a smaller value for the splitting 
5 = |((m~) — (m 2 )) (induced by a lighter spectra), and a lighter charged Higgs boson. (We notice 
that accommodating light neutral Higgs with rrih > 118 GeV is also comparatively more difficult 
in the type II SUSY seesaw.) 

Notice that in both SUSY seesaws it is fairly easy to accommodate a potential observation 
of BR(/i ei) ~ 10" 13 by MEG, taking for instance M Se esaw ~ 10 12 GeV for the type I and II 
seesaw mechanisms. 
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Figure 3: m - My 2 plane for tan/3 = 40 and A = -500 GeV, with S^ R = 0.1. On the left 
(right) panel, a type I (II) SUSY seesaw, considering degenerate heavy mediators. Contour lines 
denote values of Ar (decreasing values: positive - in association with an orange- yellow- white colour 
gradient; negative - blue gradients); solid (grey) regions are excluded due to the - requirement of 
having the correct EWSB. A green dot-dashed line corresponds to the present LHC bounds on the 
cMSSM [46]. A full green line delimits the BR(r — > ej) exclusion region, while full (dot-dashed) 
red lines correspond to the bounds on BR(i? s — > (BR(B U — > tv)). Finally, the region 

delimited by blue lines corresponds to having BR(^ — > ej) within MEG reach (current bound - 
solid line, future sensitivity - dashed line). Crosses (circles) correspond to the benchmark point 
40.1.1 (40.3.1). 



For both cases, larger values of = 0.5 can be taken, but these typically lead to conflicting 
situations with low-energy observables; lowering tan f3 can ease the existing tension, at the expense 
of also reducing Ar. We summarise this on Table [2] for simplicity in association with a type I 
SUSY seesaw. 

A few comments are in order regarding the summary of Table [2] even with a large value for 
5§±, and in the large tan fi regime, the maximum attainable values for Ar are much below the 
current experimental sensitivity, at most 2 x 10~ 6 . As mentioned before, if we further take into 
account the recent discovery of a new boson at LHC 37 with a mass around 125 GeV, and interpret 
it as the lightest neutral CP-even Higgs boson of the MSSM, the attainable values for Ar will be 
extremely small. 

In order to conclude this part of the analysis we provide a comprehensive overview of the 
constrained MSSM prospects regarding Rk, presenting in Fig. [4] a survey of the (type I seesaw) 
mSUGRA parameter space, for two different regimes of taking all present bounds (including 
the recent ones on m/j,) into account. The panels of Fig. [1] allow to recover the information that 
could be expected from the discussion following Fig. [3] for fixed values of Aq and tan /3, increasing 
^31 indeed allows to augment the SUSY contributions to Ar although, as can be seen from the 
right-panel, the constraints from BR(r — > e^f) become increasingly harder to accommodate. (Notice 
that the latter could be avoided by increasing the SUSY scale (i.e. on regions of the parameter 
space with large uiq and/or Myi) - however, and as visible from Fig. |4j in a constrained SUSY 
framework this would lead to heavier charged Higgs masses, and in turn to suppressed contributions 
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-8 
-8 
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3.1xl0~ 15 
2.9xl0~ 10 
2.0xl0~ 8 


1.09 
1.09 
1.09 


4.4 
4.4 
4.4 


1.2x10- 
1.2x10- 
3.3x10- 


-12 
-12 
-12 



Table 2: Ar and other low-energy observables for different mSUGRA points, considering a type 
I seesaw, and distinct values of S^f 1 . The values of the seesaw scale were varied from 1.3 x 10 12 
GeV to 5 x 10 10 GeV, in order to comply with the limits/future sensitivity on BR(/j — > ey). 



to Ar.) 

Although we do not display an analogous plot here, the situation is very similar for the type II 
SUSY seesaw (slightly even more constrained due to the fact that accommodating ~ 125 GeV 
is more difficult in these models |47|). 

In view of the above discussion it is clear that even taking into account all 1-loop corrections 
to the v£H + vertex, values of Ar, large enough to saturate current observation, cannot be reached 
in the framework of constrained SUSY models (and its seesaw extensions accommodating neutrino 
data). In this sense, and even though we have followed a different approach, our results follow the 
conclusions of |13| . We also stress that recent experimental bounds (both from flavour observables 
and collider searches) add even more severe constraints to the maximal possible values of Ar. 



3.2 mSUGRA inspired scenarios: inverse seesaw and L — R models 

We briefly comment here on the prospects of the inverse SUSY seesaw concerning Rk- recently, it 
was pointed out that some flavour violating observables can be enhanced by as much as two orders 
of magnitude in a model with the inverse seesaw mechanism |48| . Within such a framework, right- 
handed (s)neutrino masses can be relatively light, and as a consequence these z^j, Vr states do not 
decouple from the theory until the TeV scale, hence potentially providing important contributions 
to different low-energy processes. Nevertheless, the specific co ntribu tions to Ar are suppressed by 



,2 



a factor — f , with respect to those discussed above (see Eq. (2.14)), so that we do not expect a 
significant enhancement of SUSY 1-loop effects to Rk due to the inverse seesaw mechanism. 

For completeness (and although we do not provide specific details here), we have considered a 
specific L—R seesaw model |36| . In this framework, non- vanishing values of 5§f* can be dynamically 
generated. We have numerically verified that typically one finds 5^ < 0.01 (we do not dismiss 
that larger values might be found, although certainly requiring a considerable amount of fine- 
tuning in the input parameters). We have not done a dedicated Ar calculation for this case, but 
taking into account that the effect scales with {S^) 2 , we also expect the typical range for Ar to 
be far below the current experimental sensitivity. 
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M 1/2 (GeV) M 1/2 (GeV) 

Figure 4: mSUGRA (type I seesaw) m - Mi/ 2 plane for tan /3 = 40 and A = GeV, with 5§{ R = 
0.1 (left panel) and 5^ = 0.7 (right panel). Contour lines denote values of Ar (decreasing values: 
positive - in association with an orange-yellow- white colour gradient; negative - blue gradients). 
A full green line delimits the BR(r — > e^f) exclusion region, while full (dot-dashed) red lines 
correspond to the bounds on BR(I? S — > (BR(£? n -4 tv)). Superimposed are the regions for 

the Higgs boson mass: the dark band is for 125 < m h o < 126 (GeV) and the lighter one marks the 
region where 124 < m h o < 127 (GeV). 





m 


M 1/2 


m 2 Hi , m 2 H2 


tan f3 


°31 




(GeV) 


(GeV) 


(GeV 2 ) 






Min 





100 


-5.2 x 10 6 


40 


0.1 


Max 


1500 


1500 


-4.6 x 10 6 


40 


0.7 



Table 3: Range of NUHM parameters leading to the scan of Fig. [5] 



3.3 mSUGRA inspired scenarios: NUHM 



As can be seen from the approximate expression for Ar in Eqs. (2.14, 2.15), regimes associated 
with both large tan f3 and a light charged Higgs can greatly enhance this observable [13] (Ar oc 
tan 6 P/m A H+ ). By relaxing the mSUGRA-inspired universality conditions for the Higgs sector, as 
occurs in NUHM scenarios, one can indeed have very low masses for the H + boson at low energies. 
This regime corresponds to a narrow strip in parameter space where rrVjj ~ m/jj , in particular 
when both are close to —(2.2 TeV) 2 . In addition to favouring electroweak symmetry breaking, 
since m 2 H+ ~ | m 2 H ^ — vr? H | (even accounting for RG evolution of the parameters down to the weak 
scale), it is expected that the charged Higgs can be made very light with some fine tuning [l3] . 
In order to explore the maximal possible values of Ar, a small scan was conducted around this 
region, where changes very rapidly (see Table [3]). 

As can be verified from the left-hand panel of Fig. [5] one could in principle have semi- 
constrained regimes leading to sizable values of Rk, 0(1O~ 2 ). Once all (collider and low-energy) 
bounds have been imposed, one has at most Ar < 10 (in association with > 500 GeV). 

Moreover, it is interesting to notice that SUSY contributions to BR(B U —> tv\ which become non- 
negligible for lighter H^, have a negative interference with those of the SM, lowering the latter 
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Figure 5: Left panel: Ar as a function of the charged Higgs mass, mjj+ (in GeV). Magenta points 
have been subject to no cuts, blue points comply with the bounds on the masses (LEP+LHC), red 
points satisfy all bounds except BR(I? n —> rv) and green points satisfy all bounds. Right panel: 
BR(i? u — > tv) versus m H +. Red points satisfy only the bounds on the masses (LEP+LHC) while 
green points comply with all bounds. 



BR to values below the current experimental bound. This can be seen on the right-hand panel of 
Fig. [5] We will return to this topic in greater detail in the following subsection, when addressing 
the unconstrained MSSM. 

3.4 Unconstrained MSSM 

To conclude the numerical discussion, and to allow for a better comparison between our approach 
and those usually followed in other recent analyses (for instance |i~2j[i~4] ), we conduct a final 
study of the unconstrained, low-energy MSSM. Thus, and in what follows, we make no hypothesis 
concerning the source of lepton flavour violation, nor on the underlying mechanism of SUSY 
breaking. Massive neutrinos are introduced by hand (no assumption being made on their nature), 
and although charged interactions do violate lepton flavour, as parametrised by the C/mns matrix, 
no sizable contributions to BR(/i — > cy) should be expected, as these would be suppressed by 
the light neutrino masses. At low-energies, no constraints (other than the relevant experimental 
bounds) are imposed on the SUSY spectrum (for simplicity, we will assume a common value for 
all sfermion trilinear couplings at the low-scale, A{ = Aq). The soft-breaking slepton masses are 
allowed to be non-diagonal, so that a priori a non-negligible mixing in the slepton sector can occur. 
In order to better correlate the source of flavour violation at the origin of Ar with the different 
experimental bounds, we again allow for a single FV entry in the slepton mass matrices: 5^ ~ 0.5 
(otherwise setting all other 5^ Y = 0). 

In our scan we have varied the input parameters in the ranges collected in Table [4] We have 



also applied all relevant constraints on the low-energy observables, Eq. (3.3 3.8), as well as the 



constraints on the SUSY spectrum [5||38j. In particular we have assumed the conservative limits 

mq L R > 1000 GeV , m g > 1000 GeV . (3.13) 
Concerning the light Higgs boson mass, no constraint was explicitly imposed. We just notice here 
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Aq 
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tan/3 




other &f y 


Min 

Max 


100 
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50 
1500 


100 
2500 


1100 
2500 


-1000 
1000 


100 
2200 


100 
2500 


1200 
5000 


30 
60 


0.5 
0.5 







Table 4: Range of variation of the unconstrained MSSM parameters (dimensionful parameters in 
GeVs). Aq denotes the common value of the low-energy sfermion trilinear couplings. 



that values close to 125 GeV [37], or even larger, are easily achievable due to the heavy squark 
masses. This can be observed from the left panel of Fig. [HJ where we display the output of the above 
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Figure 6: Left panel: Ar as a function of the lightest Higgs boson mass mh (in GeV) for the range 
of parameters shown in Table [4] Red points satisfy the bounds on the spectrum (LEP+LHC), blue 
points satisfy all bounds except BR (B u — > rv) and green points satisfy all bounds. Right panel: 
mjj+ versus Aq (both in GeV), with the same colour code. Leading to both plots, the different 
input parameters were varied as in Table [4] 

scan, presenting the values of Ar versus the associated light Higgs boson mass, mh- As expected, 
no explicit correlation between mh and Ar is manifest, nor with the other (relevant) flavour-related 
low-energy bounds. For completeness, and to better illustrate the following discussion, we present 
on the right-hand panel of Fig. [6] the charged Higgs mass as a function of Aq, again under a colour 
scheme denoting the experimental bounds applied in each case. Identical to what was observed in 
Fig. [5] (notice that NUHM models correspond, at low-energies, to a subset of these general cases), 
regimes of very light charged Higgs are indeed present, in association with small to moderate 
(negative) regimes for Aq. Nevertheless, these regimes - which could potentially enhance Ar - are 
likewise excluded due to a strong conflict with BR {B u — > rv). This can be further confirmed from 
the left panel of Fig. [7j where we display the possible range of variation for Ar as a function of 
rriff+, colour-coding the different applied bounds. 

As can be seen from both panels of Fig. [7J values Ar ~ 0(1O -2 , 10 _1 ) could be obtainable, in 
agreement with Refs |llfjl4| . However, the situation is substantially altered when one takes into 
account the current experimental bounds on B decays [B u — > tu and B s — > fi + fi~) and r — > ej. 
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Figure 7: Ranges of variation of Ar in the unconstrained MSSM as a function of m#+ (left panel), 
and as a function of BR(r —> ej) (right panel). The different input parameters were varied as in 
Table |4] (notice that 5^ = 0.5). On the left panel red points satisfy the bounds on the masses 
(LEP+LHC), blue points satisfy all bounds except BR(B U — > tv) and green points comply with 
all bounds. Similar colour code on the right panel, except that blue points now comply with 
all bounds except BR {B u — > tv) and BR(r — > ej) while magenta denotes points only failing the 
bound on BR(r — > e"f). 



As is manifest from the left panel of Fig. [7j once experimental bounds - other than B u — > tv - are 
imposed, one could in principle have Ar max ~ 0(1CP 2 ); however, taking into account the limits 
from BR(B U — > tv), one is now led to Ar < 10 -3 . 

A few comments are in order regarding the impact of the different low-energy bounds from 
radiative tau decays and l?-physics observables. Firstly, let us consider the r — > ej decay: although 
directly depending on S^, its amplitude is (roughly) suppressed by the fourth power of the average 
SUSY scale, msuSY- As can be seen from Eqs. (2.14 2.15 ), Ar only depends on the charged Higgs 
mass - if the latter is assumed to be an EW scale parameter, Ar will be thus independent of msuSY 
in these unconstrained models. As such, it is possible to evade the r — > cy bound by increasing 
the soft SUSY masses, and this can indeed be seen from the right-hand panel of Fig. [7J where a 
number of "blue" points are found to lie below the BR(r — > cy) bound. 

Secondly, the B s — > decay can be a severe constraint regarding the SUSY contributions to 

Ar in the case of constrained models (see, e.g., Figs. |3]and[4|. We notice that B s — > is 
approximately proportional to Aq (see for instance |43| ) while Ar shows no such dependence: thus 
a regime of small trilinear couplings easily allows evade the B s — > fJ. + n~ bounds. 
Finally, let us discuss the B u — > tv bounds. Notice that this is a process essentially identical to the 
charged kaon decays at the origin of the Rk ratio (the only difference being that the K + meson is 
to be replaced by a B u and the e/ ji in the decay products by a kinematically allowed r), and hence 
its tree-level decay width can be inferred from Eqs. (2.1 ) and (2.3). Due to a negative interference 
between the SM and the MSSM contributions, given by the term proportional to tan 2 /3/m^ ± in 
Eq. (2.3), regimes of low lead to excessively small values of B u — > tv (below the experimental 
bound), effectively setting a lower bound for for m 2 H± (see right panel of Fig. |5j in relation to the 
discussion of NUHM models). In turn, this excludes regimes of m H + associated to sizable values 
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of Ar, as is clear from the comparison of the "blue" and "green" regions of the left panel of Fig. [7] 

In summary, we conclude that saturating the experimental bound on Rk clearly proves to be 
extremely difficult (if not impossible), even in the unconstrained MSSM, especially in view of the 
stringent constraints from B u — > tv. 

4 Conclusions 

In this work we have revisited supersymmetric contributions to Rk = T (K — > ev)/T (K — > fiv), 
considering the potential of a broad class of constrained SUSY models to saturate the current 
measurement of Rk- We based our analysis in a full computation of the one-loop corrections to the 
v£H + vertex; we have also derived (when possible) illustrative analytical approximations, which in 
addition to offering a more transparent understanding of the role of the different parameters, also 
allow to establish a bridge between our results and previous ones in the literature. Our analysis 
further revisited the Rk observable in the light of new experimental data, arising from flavour 
physics as well as from collider searches. 

We numerically evaluated the contributions to Rk arising in the context of different minimal 
supergravity inspired models which account for observed neutrino data, further considering the 
possibility of accommodating a near future observation of a [i — > ej decay. As expected from the 
(mostly) LL nature of the radiatively induced charged lepton flavour violation, type I and II seesaw 
mechanisms implemented in the cMSSM provide minimal contributions to Rk, thus implying that 
such cMSSM SUSY seesaws cannot saturate the present value for Ar. 

We then considered unified models where the flavour-conserving hypothesis on the RR slepton 
sector is relaxed by allowing a non- vanishing 5^ (e — r sector). In all models, special attention was 
given to experimental constraints, especially four observables which turn out to play a particularly 
relevant role: the recent interval for the lightest neutral Higgs boson mass provided by the CMS 
and ATLAS collaborations, BR(S S — > /i + /i~), BR(B U — > tv) and BR(r — > e^f). These last two 
exhibit a dependence on (B u — > tv) and on 5^ (r —¥ &y) similar to that of Ar. The SUSY 

contributions to Ar are thus maximised in a regime in which mjj+ and 5^ are such that the 
experimental limits for B u — > tv and r — > e"f are simultaneously saturated; in this regime one 
must then accommodate the bounds on other observables, such as irih and BR(i? s — > /x + /i _ ). For 
a minimal deviation from a pure cMSSM scenario allowing for non- vanishing values of S^, we can 
have values for Ar at most of the order of 10 -6 . In fact, the requirement of having a Higgs boson 
mass of 125-126 GeV is much more constraining on the cMSSM parameter space than, for instance 
B s — > fi^ [T~ (which is sub-dominant, and can be overcome by variations of the trilinear coupling, 
Aq). In order to have Ar ~ 0(1CU 6 ), one must significantly increase 5^ so to marginally overlap 
the regions of mh ~ 125 GeV, while still in agreement with r — > ej. 

Models where the charged Higgs mass can be significantly lowered, as is the case of NUHM 
models, allow to increase the SUSY contributions to Ar, which can be as large as 10~ 4 (larger 
values being precluded due to B u — > tv decay constraints). 

More general models, as the unconstrained MSSM realised at low-energies, offer more degrees of 
freedom, and the possibility to better accommodate/evade the different experimental constraints. 
In the unconstrained MSSM, one can find values of Ar one order of magnitude larger, Ar ~ 
0(1O~ 3 ). Again, any further augmentation is precluded due to incompatibility with the bounds 
on B u — > tv. 

However Ar ~ 0(1(T 3 ) still remains one order of magnitude shy of the current experimental 
sensitivity to Rk, and also substantially lower than some of the values previously found in the 
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literature. As such, if SUSY is indeed discovered, and unless there is significant progress in the 
experimental sensitivity to Rk, it seems unlikely that the contributions to Rk of the SUSY models 
studied here will be testable in the near future. On the other hand, any near-future measurement 
of Ar larger than 0(1CP 3 ) would unambiguously point towards a scenario different than those 
here addressed (mSUGRA-like seesaw, NUHM and the phenomenological MSSM). 

It should be kept in mind that the analysis presented here focused on the impact of LFV 
interactions. Should the discrepancy between the SM and experimental observations turn out to 
be much smaller than 10 -4 , a more detailed approach and evaluation will then be necessary. 
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A Renormalisation of the v£H + vertex 



In what follows we detail the computation leading to Eqs. ( 2.8|2,10 ), and we further refer to 41 



for a similar analysis. As expected, loop effects contribute to both kinetic and mass terms of 
charged leptons as well as to the v(.H + vertex: 



=i t L (1 + rf L ) Ul + i ?R [} + Vr ) «" h 
+ iV L {l + rf L ) 0u L - \l L {M m + rjQ £ R + h.c. 
V L (2 3 / 4 G^ /2 tan/3 M /0 + V H ) l R H + + h 



+ 



(A.l) 



Here M 10 denotes the bare charged lepton mass and the rj's correspond to loop contributions to the 
various terms. The (new) kinetic terms can be recast into a canonical form by means of unitary 
rotations of the fields (K^, K^, K£), which are then renormalised by diagonal transformations 
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Two unitary rotation matrices (R e L , R R ) are further required to diagonalise the charged lepton 
mass matrix, and one finally has 



/jold 
l R 



KUZl) 2 R{tr, 



2 r>£ anew 
n R l R > 



(A.5) 
(A.6) 
(A.7) 



In the new basis, the mass terms now read 

£ mass = _J lM 1 £r + hc 



K R Zk 



2 



R e R £ R + li.c. 



(Ai 



The above equation relates the unknown parameter M 10 with the physical mass matrix M l . Using 
the latter to rewrite the vlH + vertex one finds 



where 



C H± = v L Z H £ R H + + h.c. , 



Z H =2 3 / 4 G;/ 2 tan/Si^ (zf)~ 5 K$ K e L (z£) 5 R{ M l 
+ B$} (Zi)~ l K$ (-2^ G f tan/3^ + rf) 



(A.9) 



To one-loop order, this exact expression simplifies to 
Z H =2 3 ^G 1 J 2 tan/3 



1 + |-|) M <-^ 



+ rj 



ii 



(A.10) 



(A.ll) 
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The expressions for the 77's can be computed from the relevant Feynman diagrams (assuming zero 
external momenta): 



(4vr) 2 77 



1(4 

rn 



- (^fivlh = N^Nf^B, (o,m^) 



lap jap 1 \ ' ' . 



- (4.) = ^"S? 



L(5+)* 



( " > ' '°' m x£' m x^' m f 7 



+ ^i 2 5+) * dC, oo ( 0, 0, 0, m^o , m^± , mf_ 



D 



L(s+y 

/3a2 



C ( 0,0,0,m x o,m x ±,m~ 



~ h D fa2 ^* dC oo (0,0,0, m^o , m 2 x ± , m~ 7 



^ iv ia/3 iv ja 7 i/27/3 m X 7 



Co(o,0,0,m| 7 ,m^,m 2 o) , 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 
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with -Bo,1j Cqj Cq,0 denoting the usual loop integral functions 



B (0,x,y) 
B l (0,x,y) 



A £ + l 



x log 



yiog-^ 



x - y 
3x — y 
2(x-y) 



(A.17) 



log + 



log 



xylog-+yzlogf +zxlog 
C (0,0,0, x,y,z) = ? — ^— - — 

(a; -y){y- z) {z - x) 

dC 00 (0,0,0, x,y,z) = A e + 1 
x 2 (y - z) log 4r + y 2 (z - x) log ^ 



+ *r (x 



+ 



-y)(y - z) (z - x) 



Here d = 4 — e, /U is the regularisation parameter and A £ 
notation we followed |49|. 



(A.18) 
(A.19) 

(A.20) 

7 + log47r. For the couplings 



y) lo § 



The comparison of the above expressions with the corresponding ones derived in Ref. 41 



reveals a fair agreement; we nevertheless notice that the neutralino and chargino masses are absent 



from the analogous of Eq. (A. 12), and that the order of the arguments of B\ in Eqs. (A. 13 A. 14 



A. 15) appears reversed. Moreover, we find small discrepancies (which cannot be accounted by 



the distinct notations used) in the expressions for rf m and r]H, cf. Eq. (A. 12) and Eq. (A. 16), 
respectively. 



B SUSY seesaw models 

In its different realisations, the seesaw mechanism offers one of the most appealing explanations 
for the smallness of neutrino masses and the pattern of neutrino mixing angles. Moreover, when 
embedded in the framework of SUSY models - the so-called SUSY seesaw - the seesaw offers the 
interesting feature that flavour violation in the neutrino sector (encoded in non-diagonal neutrino 
Yukawa couplings) can radiatively induce flavour violation in the slepton sector at low-energies (8j, 
leading to potentially sizable contributions to a large array of observables. 

In what follows we briefly summarise the most relevant features of different realisations of the 
seesaw mechanism. In particular, we will consider "high-scale" seesaws, i.e., where the additional 
states are assumed to be much heavier than the electroweak scale (in association with large values 
of the corresponding couplings). 



B.l Type I SUSY seesaw 

In a type I SUSY seesaw, the MSSM superfield content is extended by three right-handed Majorana 
neutrino superfields. The lepton superpotential is thus extended as 

yyiepton = n c Y u LH 2 + E c Y l LHi + -N C M N N C , (B.l) 

where, and without loss of generality, one can work in a basis where both Y l and Mn are diagonal 
(Y l = diag(Y e ,Y^,Y T ), M N = diag (M Nl , M N2 , M Na )). The relevant slepton soft-breaking terms 



21 



are now 



VST = m l *l h + m\ Ir l* R + m\ R v R u* R + (a 1 H x h l* R 

+A U H 2 D L D R + B u D R D R + \i.c.). (B.2) 

Should this be embedded into a cMSSM, then the additional soft breaking parameters would also 
obey universality conditions at the GUT scale, (mp R )fj = m\ and (A u )ij = Ao(Y u )ij. 
In this case, the light neutrino masses are given by 

ml = -m u D T M^m v D , (B.3) 

with m v D = Y v V2 {vi being the vacuum expectation values (VEVs) of the neutral Higgs scalars, 
^1(2) = v cos (sin) /3, with v = 174 GeV), and where corresponds to the masses of the heavy 
right-handed neutrino eigenstates. The light neutrino matrix m v is diagonalized by the C/mns as 
^dmg _ jj MNS T mu u MNS ^ A convenient means of parametrising the neutrino Yukawa couplings, 
while at the same time allowing to accommodate the experimental data, is given by the Casas- 
Ibarra parametrisation |50|, which reads at the seesaw scale, M^v, 



Y v v 2 = m v D = iyjM%* R V mT g tW • (B.4) 

In the above, R is a complex orthogonal 3x3 matrix that encodes the possible mixings involving 
the right-handed neutrinos, in addition to those of the low-energy sector (i.e. C/mns) anci which 
can be parametrised in terms of three complex angles 6% (z = 1,2,3). In our analysis, we assumed 
degenerate right-handed neutrino masses and real parameters, so that the results are effectively 
independent of the choice of the Q{. 

Even under universality conditions at the GUT scale, the non-trivial flavour structure of Y v 
will induce (through the running from Mqtjt down to the seesaw scale, Mat) flavour mixing in the 
otherwise approximately flavour conserving soft-SUSY breaking terms. In particular, there will 
be radiatively induced flavour mixing in the slepton mass matrices, manifest in the LL and LR 
blocks of the 6x6 slepton mass matrix; an analytical estimation using the leading order (LLog) 
approximation leads to the following corrections to the slepton mass terms: 

(Am|) y = -^(3ml + Al)(Y^LY^, 

(Ami),. ~ ; L M ee log (^~) hi ■ (B.5) 



The amount of flavour violation is encoded in the matrix elements (Y u 'LY v )ij of Eq. (B.5) 



B.2 Type II SUSY seesaw 

The implementation of a type II SUSY seesaw model requires the addition of at least two SU(2) 
triplet superfields [51] . Should one aim at preserving gauge coupling unification, then complete 
SU(5) multiplets must be added to the MSSM content. Under the SM gauge group, the 15 
decomposes as 15 = S + T + Z, where S ~ (6, 1, —2/3), T ~ (1, 3, 1) and Z ~ (3, 2, 1/6). In the 
SU(5) broken phase (below the GUT scale), the superpotential contains the following terms: 
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Wii 



V2 



Yt LTi L + Ys D S D) + Y z D Z L 



+ Y d D c QH l + Y u U C QH 2 + Y l E c L U x 
+ ~ (Aiffi fi #1 + A 2 tf 2 T 2 tf 2 ) + M T fi f 2 
+ M z Zi Z 2 + M s 5i S 2 + HH1H2, 



(B.6) 



where we have omitted flavour indices for simplicity (for shortness we will not detail the soft 
breaking Lagrangian here, see e.g. |5l]). After having integrated out the heavy fields, the effective 
neutrino mass matrix then reads 

TT Vo A 2 



■m, 



-Y T . 



(B.7) 



2 M T 

As occurs in the type I seesaw, LFV entries in the charged slepton mass matrix are radiatively 



induced, and are proportional to the combination YiYp 51 ; for example, the LL block reads 



(A m !)« oc (Y^ Y T )ij ~ ( 



A 2 v\ 



13 



fBi 
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